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ABSTRACT 

We present extensive calculations on the efficiency of grain alignment by the 
Davis-Greenstein mechanism. We model the grains as oblate spheroids with 
arbitrary axis ratios. Our description of the grain dynamics includes (i) magnetic 
dissipation and the inverse process driven by thermal fluctuations in the grain 
magnetization; (ii) gas-grain collisions and thermal evaporation of molecules from 
the grain surface; (iii) the transformation of rotational energy into heat by the 
Barnett effect and the inverse process driven by thermal fluctuations; and (iv) 
rapid Larmor precession of the grain angular momentum about the interstellar 
magnetic field. For ordinary paramagnetic grains, we calculate the Rayleigh 
reduction factor, R, for >1000 combinations of the 3 dimensionless parameters 
which characterise the alignment. For superparamagnetic grains, we calculate R 
from an exact analytic solution for the relevant distribution function. Our results 
are compared with classical studies of DG alignment, which did not include the 
Barnett effect. We calibrate the accuracy of a recently-proposed perturbative 
approximation, which includes the Barnett effect, and show that it yields R 
values with a mean error of ~17%. 

Subject headings: dust - - polarization — ISM: continuum — ISM: magnetic 
fields 



1. Introduction 

Nearly half a century after the discovery of interstellar polarization, the alignment of 
interstellar dust grains remains poorly understood. One candidate for the mysterious align- 
ment mechanism is a process described by Davis and Greenstein (1951, henceforth DG51), in 
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which paramagnetic dissipation aligns the grain angular momentum, J, with the interstellar 
magnetic field, B. The alignment of the grain axes with B requires a correlation between 
J and the axial directions. In the original formulation of DG51, the required correlation 
is produced by the torques due to paramagnetic dissipation and gas damping. The Davis- 
Greenstein (DG) mechanism was subsequently modified to describe ferromagnetic grains 
(Spitzer & Tukey 1951) and to include the important disorienting effects of thermal fluctu- 
ations in the magnetization (Jones & Spitzer 1967, henceforth JS67). The "classical age" of 
grain alignment research culminated with extensive quantitive investigations based on the 
Fokker-Planck equation (JS67) and Monte Carlo methods (Purcell 1969; Purcell & Spitzer 
1971, henceforth PS71). These studies concluded that the DG mechanism can explain the 
polarization observed toward diffuse clouds if the grains are superparamagnetic, but requires 
magnetic fields > 10~ 5 G if the grains are composed of ordinary paramagnetic substances. 

The classical era was followed by a twenty-year period with little research on the DG 
effect, during which interest focussed on an alternative mechanism proposed by Purcell 
(1979, henceforth P79). In Purcell's mechanism, magnetic dissipation aligns the angular 
momenta with B, as in DG51, but rotation at superthermal energies stabilizes the grains 
against disorientation by the gas particles. P79 also showed that the "internal alignment" 
of the grain axes with J is determined by the conversion of rotational energy into heat via 
frictional processes internal to the grain (e.g., the Barnett effect). Internal dissipation drives 
an isolated grain toward the state of minimum energy consistent with angular momentum 
conservation. For rotational energies much larger than the equipartition energy, the result 
is almost perfect alignment between J and the principal axis of largest rotational inertia. 

Research on the DG mechanism has been reinvigorated by studies of polarization toward 
molecular clouds, where the physical conditions appear unfavorable for Purcell's mechanism 
(and other mechanisms, see Hildebrand 1988, 1996 and references therein). Roberge, De- 
Graff and Flaherty (1993, henceforth RDGF93) modified the DG effect to include internal 
dissipation in the regime where the internal alignment is perfect. Lazarian (1995, henceforth 
L95) presented an accurate analytic solution for the alignment efficiency in this regime. De- 
tailed studies of internal dissipation (Lazarian 1994; Lazarian & Roberge 1997, henceforth 
L94 and LR97, respectively) showed that internal alignment becomes imperfect at thermal 
rotation energies and produced an accurate quantitative theory of the dissipation mecha- 
nism. The modification of the DG effect to include imperfect internal alignment was carried 
out by Lazarian (1997, henceforth L97), who gave an approximate analytic solution of the 
Fokker-Planck equation. 

The approach of L97 replaces certain random variables in the Fokker-Planck equation 
with their mean values. Though physically well motivated, this is an ad hoc approximation 
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whose accuracy can only be determined by comparison with calculations of known accuracy. 
In this paper we present such calculations. Our plan is as follows: The dynamical processes 
and modelling assumptions included in our calculations are summarised in §2. In §3 we 
describe an adiabatic approximation (Roberge 1997, henceforth R97) that yields accurate 
solutions of the Fokker-Planck equation along with a rigorous procedure for estimating the 
errors. In §4 we evaluate these solutions numerically for ordinary paramagnetic grains (§4.1) 
and analytically for superparamagnetic grains (§4.2). We compare our results with the 
classical studies of DG alignment in §5.1 and with the results of L97 in §5.2. The status of 
the DG mechanism, including outstanding problems, is discussed in §6 and our results are 
summarised in §7. 



We model the grains as oblate spheroids with semiaxes a parallel to the symmetry axis 



The components of the inertia tensor for rotation parallel and perpendicular to z will be 
denoted I\\ and I±, respectively. Our model grains are composed of an unspecified solid 
with uniform density p s , temperature T s , and frequency-dependent magnetic susceptibility 
X (u>) = x' i 00 ) + l x" They are located in a static gas of particles with mass m, number 
density n, and kinetic temperature T g , plus a uniform magnetic field, B. 



We are interested ultimately in the distribution of the grain orientations in an inertial 
frame fixed to the gas. Let {x, y, z} be the basis of the inertial frame, with z parallel to 
B and the other basis vectors oriented arbitrarily in the plane perpendicular to B (Fig. 2). 
The orientation of a in the inertial frame changes due to (i) the motion of the grain's axes 
with respect to J caused by rotation and the Barnett effect, plus (ii) the motion of J in 
the inertial frame caused by Larmor precession and the DG effect. It is therefore natural to 
specify the orientation in terms of the four angles defined in Figures 1 and 2, so that (9, i/j) 
gives the orientation of a with respect to J and (f3, 0) gives the orientation of J with respect 



2. Modelling assumptions 



2.1. 



Grain model 




2.2. Dynamical processes 



to B. 



The dynamical timescales for these variables are listed in Table 1. The rotational motion 
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of a spheroid causes a to precess about J with period ~ t rot (for reasonable grain shapes) 
and the interaction of the grain magnetic moment with the interstellar magnetic field causes 
J to precess about B with period t^. The smallness of the rotation and Larmor periods 
compared to all other timescales of interest insures that and 4> are uniformly distributed. 
Thus, only the joint distribution of 9 and (5 needs to be calculated in practice. 

We assume that internal dissipation is dominated by the Barnett effect]], which there- 
fore determines the the motion of $.0 The systematic and random motions of 9 produced 
by Barnett dissipation and thermal fluctuations in the Barnett magnetization (collectively 
termed "Barnett relaxation") have dynamical timescale 

^ 2/ ll 

tBar = VKh\h-\)P (2_1) 

(RDGF93), where V is the grain volume, h = I\\/I±, K = x" /uo and /i is the magnetogyric 
ratio of the orientable spins or orbits responsible for the Barnett magnetization. The Barnett 
time depends on the angular momentum; the value given in Table 1 was computed for a grain 
with J equal to 



Jth = yJI\\kT s , (2-2) 

a typical value for thermal rotation. 

We assume that the motion of (3 is caused by (i) the systematic and random torques 
produced by gas-grain collisions and thermal evaporation of molecules from the grain surface 
("gas damping"); and (ii) the systematic and random torques produced by paramagnetic or 
superparamagnetic relaxation ("magnetic damping"). The gas damping time is the timescale 
for J to be reduced by rotational friction with the gas. For a spheroid rotating about its 
symmetry axis, 



^gas — A i~2 t_a_. -n ^) 



3/|| 

(RDGF93), where = ^2kT g /m is the gas thermal speed and Ty is a dimensionless 
coefficient that depends weakly on the grain shape [see eq. (|A7|) 1- The magnetic damping 



1 A recent study by Lazarian & Efroimsky (1999) shows that, for grains with large axis ratios and grains 
produced by agglomeration, internal dissipation is dominated by inelasticity. However, this has no effect on 
the conclusions of the present study, which only requires the timescale for internal dissipation to be much 
smaller than the timescale for external interactions. 

2 In principle, we should also include gas damping: the impulsive change in J due to a gas-grain collision 
or evaporation from the grain surface changes 9. In omitting these effects we introduce errors in the 9 
distribution that, for typical molecular clouds, are of order tear/igas ~ 10~ 3 . 
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time for a spheroid rotating about its symmetry axis is 

*mag = -J^ffl- ( 2 - 4 ) 

The choice of a particular rotation axis for the purpose of defining characteristic timescales 
obviously involves no loss of generality, but must be taken into account when comparing 
results from different papers. In particular, PS71 defined t gas and t mag with respect to 
rotation about the transverse axis. To compare our results with those of PS71, it is only 
necessary to note that 5 m = (T±/T\^j 5, where T± is a dimensionless function of the grain 
shape [see eq. (|A8|)], 5 m is our magnetic damping parameter [see eq. (|3-17|) 1 and 5 is the 
magnetic damping parameter of PS71. 



2.3. Characterisation of the alignment 

Our objective is to calculate the quantity which characterises the mean axial align- 
ment for an ensemble of identical grains. Because different definitions of this quantity have 
appeared in the literature, due to differing assumptions about the grain dynamics, a brief 
review of this subject is appropriate. 

Consider an electromagnetic wave propagating along the z° axis of the Cartesian "ob- 
server frame" defined in Figure 3. The transfer equations for the Stokes parameters depend 
on the cross sections^ C° and C° for linearly polarized waves with the electric vector, E, 
along the x° and y° directions (e.g., see Martin 1974, Lee & Draine 1985). To calculate these 
"observer frame" cross sections, one transforms the components of E to a frame aligned with 
the principal axes of the grain and takes the appropriately-weighted sum of the cross sec- 
tions, Cii and C±, for E polarized along the grain axes. When the transformation is carried 
out and the resulting expressions are averaged over ip and 0, one finds that the mean cross 
sections for an ensemble of oblate spheroids are 

C° = C avg + 1r(c ± - C||) (l - 3 cos 2 C) (2-5) 

and 

C° y = C &vg + -R (C ± -C,|), (2-6) 

where C avs = (2C ± + C||)/3 is the effective cross section for randomly-oriented grains. Grain 
alignment is characterised by the Rayleigh reduction factor, 

R=(G (cos 2 6) G (cos 2 (3) ) , (2-7) 



3 We use the term "cross section" generically to refer to the extinction and phase advance cross sections. 
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where 9 is the angle between the axis of largest rotational inertia (i.e., z b for a homogeneous 
oblate spheroid), (3 is the angle between J and B, 

G(x) = -(x~-) (2-8) 



2 V 3 

and () denotes the ensemble average. Thus, our objective is to calculate R. 

Two related "alignment measures" which frequently appear in the literature are 

Qx = (G (cos 2 e)) = I ((cos 2 6) - I) , (2-9) 



2\\ I 3 
which measures the mean alignment of a with J, and 

Q j = (G (cos 2 0)) = | ((cos 2 p) - ~) , (2-10) 

which describes the alignment of J with B. The Rayleigh reduction factor cannot be de- 
termined uniquely from Qj and Qx without making additional approximations because in 
general 9 and f3 are correlated. For example, if one assumes that Barnett dissipation aligns 
a perfectly with J (e.g., see Lee & Draine 1985; RDGF93; L95), then G{6) = 1 and R = Qj. 
This assumption is appropriate when the grains rotate with energies 3> kT s (P79) but a poor 
approximation for thermal rotation, the case of interest here (L94; LR97). We make no a 
priori assumptions about a possible relationship between R, Qx, and Qj. 



Mathematical formulation 



3.1. The Fokker-Planck equation 



Although we are interested only in the distribution of 9 and j3, it is easier for practical 
reasons to calculate the joint distribution of 9 and all 3 angular momentum coordinates. 
Let (J x , J y , J z ) be the Cartesian coordinates of J in the inertial frame. Define the distribu- 
tion function, / (9, J x , J y , J z ), so that / (9, J x , J y , J z ) d9 dJ x dJ y dJ z is the joint probability of 
finding a oriented in the infinitesimal interval d9 centered at 9 and the angular momentum 
in dJ x dJ y dJ z centered at (J x , J y , J z ). Then / satisfies the Fokker-Planck equation, 



df 
dt 



d_ 

09 



A e f - ^ (Beef) 



d 
dJ k 



_d_ 



(Bkkf) 



(3-1) 



where repeated indices are summed over x, y, z and we have anticipated the result [see 
eqs. ( p-19|) - (|3-21] )1 that B is a diagonal tensor. The diffusion coefficients Ae and Bgg, which 
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describe the variations in due to Barnett relaxation, were derived in LR97. However they 
are not needed in the approximation described below (see §3.2). The coefficients and 
describe the analogous changes in the angular momentum components due to the combined 
effects of gas- and magnetic damping. They are given for reference in Appendix A. The 
angular momentum diffusion coefficients are functions of 0, {Jk} and physical parameters 
such as the grain temperature and magnetic field strength. We assume that all physical 
parameters are constant over times much longer than the timescales in Table 1, so that only 
the steady solutions of eq. (|3-1|) are of interest. 



3.2. Adiabatic elimination of 6 

According to Table 1, the dynamical timescale for is typically ~3 orders of magnitude 
shorter than the timescales for the angular momentum components.^ It is possible to exploit 
this disparity to greatly simplify the calculation of /. In particular, R97 showed that 

/ (0, J x , J y , J z ) = f- mt (0 | J) /act (-4, Jy, J x ) + 0(e), (3-2) 

where / int {0 \ J ) can be calculated from equilibrium thermodynamics,^ / ext satisfies a simpli- 
fied Fokker-Planck equation described below and e ~ is the ratio of tear to the shortest 
timescale for the motion of J. 

LR97 showed that the conditional distribution is 

/int (0 I J ) = C e sin exp (~e sin 2 0) , (3-3) 
where Cg is a normalization constant and 

(ft -i)j» 

2hkT s K ' 



Expression ( [3-3j ) is the thermal equilibrium distribution for an isolated grain with fixed 
angular momentum. It depends parametrically on J as the notation /i nt (9\J) suggests. This 
result has an obvious physical interpretation: as the angular momentum slowly changes, the 
rapid evolution of allows f- int to relax adiabatically to the equilibrium distribution that 
would obtain if J were frozen at its instantaneous value. Consequently, we refer to expression 
as the adiabatic approximation. 



4 There is no need to modify this statement if the grains are superparamagnetic, since superparamagnetism 
would reduce i mag and tBar by the same factor. 

It follows that /i n t is independent, to 0(e), of the particulars of the relaxation mechanism. 
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The statistics of f\ nt are functions of J 2 . Examples that appear frequently in the fol- 
lowing discussion are the conditional mean values of sin 2 9, 



a 



r sin 2 9f int {9 \J)d9, 
Jo 



cos 2 9, 



and the related alignment measure 



k = 1 — er, 



qx = 



K — 



(3-5) 
(3-6) 

(3-7) 



They are plotted in Figure 4. The alignment of a with J is different for each grain in 
the ensemble, with qx increasing from zero for grains with J = to unity for grains with 
J 2 ^> (h — l)L\kT s . The mean alignment for the ensemble is described by Qx, the average 
of qx over the angular momentum distribution. 

In the adiabatic approximation, the angular momentum distribution satisfies the Fokker- 
Planck equation 

d r v . i Id 2 



dJ k 



A k f, 



ext 



+ 



2d J 2 



Bkk f, 



ext 



0. 



(3- 



where the coefficients and are obtained by averaging the ^-dependent coefficients in 
expression (|3-1|) in the obvious way Thus, 



Ai 



and 



B 



kk 



d9A k f- mt (9\J) 



d9 B kk f int (9\J). 



(3-9) 



(3-10) 



With 9 thus eliminated, the grain alignment problem reduces to solving equation ( |3-8; ) for the 
angular momentum distribution. Once this is accomplished, the Rayleigh reduction factor 
can be found by evaluating 



R = J d 3 J f cxt (J x , J y , J z ) q x (J 2 ) G (cos 2 (3) , 
where J denotes the definite integral over all of angular momentum space. 



(3-11) 



3.3. Dimensionless parameters 



We minimize the number of independent parameters in the problem by writing A and B 
in terms of dimensionless units, with angular momentum measured in units of J t h arid time 
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measured in units of t gas . Henceforth it is implied that all quantities are in dimensionless 
units unless it is explicitly stated otherwise. 

After evaluating expression ( |3-9|) and writing the result in dimensionless form, one finds 
that the mean torque has components 



and 

The damping rates, 
and 



A - —h T 

-^x n xy °xi 


(3 


A — — h 1 

V *^xy °y 


(3 


A z k z J z . 


(3 



k xy = 1 + K + [A + (h - l)6 m ] a (3-15) 



k z = 1 + Xa, (3-16) 
are functions of J 2 (through a). They depend also on the magnetic damping parameter, 

5 m = (3-17) 

''mag 

and the dimensionless shape parameters h — 1 and 

A = hT ± /T ll - 1 (3-18) 

(see Fig. 5). 

The evaluation of the diffusion tensor yields 

B xx = (1 + T s /T g ) (1 - 1Vxy ) + 2 (T s /T g ) 5 m , (3-19) 

Byy = (3-20) 

and 

B zz = (1 + T s /T g ) (1 - 1Vz ) + 2 (T s /T g ) S m , (3-21) 

where 

7 = i- r ± /r,| (3-22) 

is a dimensionless function of the grain shape (Fig. 5). The diffusion tensor depends on J 2 
and cos 2 (3 via the coefficients 



T) xy (j 2 ,COS 2 /?) EE i 



l + cos 2 /5+ -all -3cos 2 /? 



(3-23) 
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and 

rj z ( J 2 , cos 2 (3) = 1 - cos 2 /3 - -a (l - 3 cos 2 p) . (3-24) 

We conclude that the diffusion coefficients, and therefore R, depend on T s /T g , 5 m and 3 
dimensionless functions of the grain shape. We will parameterize our solutions by T s /T g , 5 m 
and a/b. 



4. Results 

4.1. Ordinary paramagnetic grains 

Equation (|3-8|) must be solved numerically in general. We do this using an algorithm 
described in RDGF93, to which we refer the reader for a thorough description of the method. 
In brief, we exploit the mathematical equivalence of eq. ( |3-8| ) and the coupled Langevin 
equations 

dJk = Akdt + y/BkkdWk (no summation on k), (4-1) 

where {(iW 7 *;} are independent Gaussian random variables with variance dt. We integrate 
these equations numerically as an initial value problem to generate simulations of the time- 
dependent angular momentum components and compute statistics of the angular momentum 
distribution, such as R, by equating time- and ensemble averages. Theorems insure that our 
statistics converge to their exact values in the limits At — > and T — > oo, where At is 
the time step and T is the total averaging time. The benchmark tests and "production" 
calculations reported in this paper use At = 10~ 3 and T = 10 5 dimensionless units unless it 
is stated otherwise. 

Because we are interested only in statistics of the steady-state distribution, the initial 
conditions on eq. (|4-1|) are arbitrary. Our results were obtained by initially setting the 
magnitude of J to unity and choosing its direction randomly from an isotropic distribution. 
To allow our simulated grains to "forget" these initial conditions, we integrated the Langevin 
equations for 1000 dimensionless time units before commencing our calculation of the time 
averages. 

Useful checks on our numerics are provided by special cases where analytic solutions 
for certain statistics exist. For example, Qj can be calculated exactly for the special case 
a/b = 1 (i.e., for spherical grains; see JS67, PS71). Let Qj (a random variable) denote the 
estimate of Qj predicted by one numerical integration with T = 10 5 and SQj = Qj — Qj 
be the corresponding error. In Figure 6 we plot the error distribution obtained from 100 
numerical integrations with a/b = 1 and the other parameters set arbitrarily to 5 m — 1 and 
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T s /T g = 0.5. The mean and rms errors are 5Qj = —1.5 x 10~ 4 and (SQ 2 ^ = 1.2 x 10~ 3 , 
respectively. We infer that the uncertainty in Qj is dominated by statistical fluctuations in 
the simulations rather than systematic errors in the numerics. 

In the limit T s /T g — > 0, the alignment of a with J becomes perfect and R — > Qj. L95 de- 
veloped a perturbative approximation which calculates Qj in this regime. Figure 7 compares 
results obtained with our numerical method for the case T s /T g = 0, a/b = 2/3 (symbols) and 
results obtained with the perturbative approximation in fifth order (solid curve). The mean 
and rms discrepancies between the 49 numerical points and the corresponding predictions of 
perturbation theory are 5Qj = 1.2 x 10~ 2 and (5Q 2 j^ = 1.3 x 10~ 2 . These discrepancies 
should be interpreted as upper limits on the errors in the numerical results, since they are due 
to the combined effects of errors in the numerics and the finite accuracy of the perturbation 
method. 

Another check on the numerics is provided by the case T s /T g = 1, where the angular 
momentum distribution is Maxwellian and Qx can be calculated exactly (see JS67). Figure 8 
shows the errors in Qx for 100 trials with T s /T g = 1 and the other parameters arbitrarily 
chosen to be a/b = 0.5 and 5 m = 1. The mean and rms errors are 5Qx = —5 x 10~ 4 and 

[oQx) = 5 x 10 , respectively. 

There do not appear to be any benchmarks that would permit us to calibrate the errors 
in R — the quantity of ultimate interest — for more general cases. However rough estimates 
can be obtained by assuming (albeit somewhat unrealistically, see the discussion below) that 
6 and f3 are uncorrelated. With this assumption, one can easily show that the mean square 
error in R satisfies 

5R? = Q 2 x 6Q^ + Q\ x 5Q 2 , + terms, (4-2) 

where the omitted terms involve higher-order products of the small quantities 5Qj and 5Qx- 
If we use the fact that Qj and Qx cannot exceed unity and assume, reasonably, that the 
values of 5Q 2 and 8Q 2 X from the benchmarks above are typical, then 

(lR2) 1/2 £ 1 x 10^ 3 . (4-3) 

Thus, we estimate a typical error in R to be a few units in the third decimal place. 

We have computed the alignment of ordinary paramagnetic grains for 1259 parameter 
combinations on the uniform grids a/b = 0.1(0.2)0.9, T s /T g = 0(0.1)0.9, and log<5 m = 
— 1 (0.1) 1.5. Our results are given in Tables 2-6.| We have deleted the table entries whenever 



These data are also available in electronic form; inquiries should be directed to WR. 
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R < 5 x 10~ 3 , consistent with our error estimate. Davis- Greenstein alignment with R values 
in the omitted range may be relevant to polarization in the 2175 A extinction feature, which 
has been attributed to magnetically-aligned graphite grains with R = 1-2 x 10~ 3 (Wolff 
et al. 1997). This hypothesis is examined in a forthcoming paper (Roberge & Karcz, in 
preparation), where DG alignment of graphite grains is studied using perturbation theory. 

The dependence of various alignment measures on the parameters is shown in Figures 9- 
11. The effects of Barnett relaxation on the grain dynamics are illustrated by the graphs of 
Qx- Barnett dissipation generally enhances the alignment of a with J, in the sense that Qx 
is typically several times larger than the Qx value for a Maxwellian distribution. Thermal 
fluctuations in the Barnett magnetization account for the systematic dependence of Qx on 
the parameters. These fluctuations cause random changes ~ kT s in the rotational energy, 



2 



F ^ 

rot ~ m 



1 + (/i - 1) sin 2 . (4-4) 



The resulting variations in 9 inhibit the alignment of a with J unless 

2hkT< 



(h-l)J 2 



« 1, (4-5) 



that is, unless kT s is much smaller than the difference between the minimum (6 = 0) and 
maximum (0 = 7r/2) rotational energies. It follows from equation ( J4-5|) that Qx should 
decrease as the grains become more spherical (so that h — > 1), as the magnetic damping 
parameter increases (so that the mean value of J 2 decreases) and as the dust temperature 
increases. These trends are present in Figures 9, 10 and 11, respectively. The graphs of Qx 
also provide additional checks on our calculations. Thus, Qx vanishes in the limit a/b — > 1 
(spherical grains disoriented completely by an infinitesimal fluctuation) and Qx approaches 
the Maxwellian solution as T s /T g — > 1. 

The graphs of Qj provide similar insight. The alignment of J with B is nearly indepen- 
dent of the grain shape, Qj being very close to the "spherical solution" for all values of the 
axis ratio (Fig. 9), magnetic damping parameter (Fig. 10) and dust-to-gas temperature ratio 
(Fig. 11). This is clearly a consequence of the weak dependence of the angular momentum 
diffusion coefficients on the grain shape and orientation.[] The monotonic increase of Qj 
with 5 m (Fig. 10) reflects the increasing efficiency of magnetic damping as the ratio t gas /t mag 
increases. Of course, the alignment should saturate, with Qj approaching some asymptotic 
limit as 5 m — > oo. However, the largest 5 m values in Fig. 10 are too small to see this limiting 



This weak dependence made the perturbative approach of L95 and L97 possible. 



13 



behavior, which is manifested only for the very large S m values characteristic of superparam- 
agnetic grains (see §4.2). The rapid decline of Qj with increasing dust temperature (Fig. 11) 
illustrates the effects of thermal fluctuations in the magnetization, whose interaction with 
B produces random torques that tend to disorient J (JS67). As T s /T g — > 1, a balance 
is reached between the dissipation of rotational energy into heat by magnetic dissipation 
and the transfer of energy to rotation by the fluctuations. In this limit, the distribution of 
the transverse angular momentum components becomes Maxwellian and the alignment of J 
vanishes. Our solutions also satisfy this benchmark. 

Our numerical code also calculates the correlations between 9 and (3. Define the dimen- 
sionless "correlation function" by 

= (cos 2 9 cos 2 g) 
9 ~ (cos 2 9) (cos 2 9) ' 1 ' 

In Figure 12 we plot p vs. 5 m for the case a/b = 2/3 and T s /T g = 0.2. The correlations 
appear to be a minor effect, in the sense that (cos 2 9 cos 2 (3) never differs by more than a 
few percent from the product (cos 2 9) x (cos 2 9) . However the apparent smallness of p is 
misleading. Let 

R = QxQj (1 + r) , (4-7) 

so that t = corresponds to the situation where the correlations vanish. It is easy to show 
from expressions (|2-7|) - (|2-10|) that p and r are related by 

1 + 2Q X + 2Qj + AQ x Qj 



4QxQ„ 



P, (4- 



so that r ^> p if QxQj *C 1. For example, consider the case 5 m = 1, T s /T g = 0.2 and 
a/b = 2/3. The alignment measures are Q x = 0.1573, Qj = 8.761 x 10~ 2 , R = 2.376 x 10~ 2 
and the correlation function is p = 0.025. Thus, R is a factor pa 1.7 larger than the product 
QxQj even though p 1. By this measure the correlations are large. In fact, this example 
is typical. For the 239 data points in Figures 9-11, the ratio R/QxQj varies from 1.3 to 5.0 
with a mean value of 1.7. In contrast, p, never exceeds 0.1. We conclude that the Rayleigh 
reduction factor is very sensitive to small correlations between 9 and j3. 

Because R is sensitive to the correlations, it is important to ask whether they are 
real or merely an artifact of the numerics. The benchmark calculations in Figs. 6-8 do 
not address this question because they describe special cases where the correlations vanish. 
Nevertheless, it is easy to see that the behavior of p in Fig. 12 is qualitatively correct: In 
the adiabatic approximation, the correlation function is merely a statistic of the angular 
momentum distribution, / ext . We cannot check the S m dependence of / ext directly because 
our numerical integration scheme yields only statistics of the various distributions. However, 
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there can be no qualitative differences between the 5 m dependence of f ext for grains with a/b = 
2/3 and / ext for spherical grains (which can be calculated exactly, see PS71). The dashed 
curve in Figure 12 is an approximation to the correlation function obtained by replacing f ext 
with the exact solution for spheres with 5 m — 1 and T s /T g = 0.2.Q The qualitative similarity 
of the two curves in Fig. 12 shows that the correlations are not a numerical artifact. An 
additional check on the correlations is discussed in §4.2, where we show that numerically- 
computed R values converge to an analytic solution for R that becomes exact in the limit 
8 m — > oo. Given the sensitivity of R to the correlations, the convergence can only occur if 
the numerical treatment of the correlations is highly accurate. 



4.2. Superparamagnetic grains 



The magnetic susceptibility of interstellar grains can be enhanced to "superparamag- 
netic" values if the grains contain small clusters of ferromagnetic or ferrimagnetic atoms 
(JS67; Draine 1996; Draine & Lazarian 1998, henceforth DL98). The existence of super- 
paramagnetic grains is consistent with the wavelength dependence of continuum polarization 
(Mathis 1986) and is supported by the presence of small FeNi and FeNiS inclusions in the 
GEMS component of interplanetary dust^ (Bradley 1994; Martin 1995; Goodman & Whittet 
1995). The alignment of superparamagnetic grains is therefore of considerable interest. 

Superparamagnetism corresponds to the case 1 <C 5 m ^ 10 5 , a regime where numerical 
calculations of the alignment are impractical. The numerical difficulties are caused by the 
disparity between the timescales for gas damping (=1 in dimensionless units) and magnetic 
damping (=5~ 1 ). (That is, the Langevin equation becomes "stiff.") To produce accurate 
statistics, numerical integrations must be carried out with a timestep much smaller than the 
smallest dynamical timescale and an averaging time much larger than the largest timescale. 
It follows that the number of arithmetic operations per R value increases in proportion to 
5 m and a numerical investigation of superparamagnetic alignment would require orders of 
magnitude more computation than the present study. 

However, the disparity between t gas and t mag can be turned to advantage using the 
adiabatic elimination technique. In Appendix B we carry out the analysis to show that 

/ext = Uy (4, Jy\J z ) f Z (J Z ) + O (C) ■ (^) 



8 That is, we computed (cos 2 6), etc using expressions (3-5)— (3-7) for nonspherical grains with a/b = 2/3 
and T s /T g — 0.2 and averaged these expression over the angular momentum distribution for spheres. 

9 DL98 claim, however, that not more than 5% of interstellar iron can be in the form of metallic inclusions. 
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The angular momentum distribution functions, f xy and f z , have closed-form solutions. For 
the transverse components, we find 

fxy (J x , Jy\Jz) = C xy exp (-$) , (4-10) 
where C xy (a function of J z ) is determined by the normalization condition 

/+oo r+oo 

/ f xy {J x ,Jy\J z )dJ x dJy = l. (4-11) 
-co J — oo 

The argument of the exponential is 

$ (J x , J, | J z ) = (^j [l + (h-l)a (J')] df ± , (4-12) 

where J± = (j% + J y ^j ^ . The distribution of J z is 

/ 2 (J 2 ) = C 2J Bexp(-^), (4-13) 
where C z is a normalization constant determined by 

/+oo 

f z (J z )dJ z = l (4-14) 
-oo 



and 



f 


" -2i 2 " 


/o 





rfj^. (4-15) 



The coefficients A and -B are functions of J 2 determined by averaging the diffusion coefficients 
in equation (|3-8|) over the transverse angular momentum components, thus: 



A, = ££££A M f xy dJ x dJ y (4-16) 

and 

B zz = ££f±™B„ f xy dJ x dJ y . (4-17) 

This completes the solution. 



One can easily verify that solution (|4-9|) reduces to the exact analytic solution for spher- 
ical grains (PS71) in the appropriate limit. An additional check is provided by comparing R 
values obtained analytically from eq. (4-9) with numerical solutions of the Langevin equa- 



tions for large 8 m . The results of a few tests with 5 m = 100 are given in Table 7. The 
numerical R values were obtained by integrating the Langevin equations with At = 10~ 5 
and T = 10 5 dimensionless units.Q According to expression ( |4-9| ), the numerical and analytic 



10 



It took 3.5 CPU days on a Sun 4 to compute each of these points. 
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results should agree to within ~1% if the discrepancies are due mainly to the accuracy of the 
adiabatic approximation rather than errors in the numerics. The discrepancies are 0.3-2%, 
consistent with this interpretation. The fact that the numerical R values are systematically 
smaller than the analytic results reinforces this conclusion. 

The Rayleigh reduction factor for superparamagnetic grains depends only on a/6 and 
T s /T g . Values of R can be calculated straightforwardly, in principle, from expressions ([4-91 ) 
and (|3-11 ). However this involves the numerical evaluation of certain multiple integrals and 



a large amount of computation. We have therefore tabulated the R values for superparam- 
agnetic grains in Table 8. 



5. Comparison with earlier work 

5.1. Classical studies 

The objective of this paper has been to provide accurate calculations on the Davis- 
Greenstein effect incorporating all physical processes now believed to be relevant. The im- 
portant issue of how our calculations bear on the classical "grain alignment problem" is 
deferred to a future paper, where we make a comprehensive comparison between the obser- 
vational constraints, our calculations on the DG effect and the calculations from forthcoming 
papers on other alignment mechanisms. Nevertheless, a brief comparison of the our results 
with previous studies of DG alignment is appropriate. 

Until recently, the standard work on DG alignment was the Monte Carlo study of PS71, 
who described the alignment of oblate and prolate grains including gas-grain collisions, evap- 
oration of molecules from the grain surface, magnetic dissipation, and thermal fluctuations 
in the grain magnetization. PS71 did not include Barnett relaxation, whose relevance to 
grain alignment was not appreciated until the work of P79. The nature of the Monte Carlo 
method used by PS71 forced them to adopt some computational expedients^ but otherwise 
their simulations model the alignment with complete verisimilitude. Thus, the essential dif- 
ference between our calculations and those of P71 is the inclusion of Barnett relaxation in 
the present study. 

In Figure 13 we compare values of Qx, Qj and R obtained with our Langevin code (solid 
curves) and values for corresponding cases tabulated in PS71 (symbols). Not surprisingly, 



For example, they modelled the grain surface with just 4 patches, scaled up the gas/grain mass ratio from 
a realistic value of ~ 10 -12 to "something like 10~ 2 " and represented the Maxwellian velocity distribution 
of the gas particles by two delta functions. 
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there is good agreement on Qj, the results from both studies being very close to the "spherical 
solution" for reasons discussed in §4.1. However there are large differences in the solutions 
for Q X - In PS71, the motion of a in the body frame is caused by gas damping. Thus, the 
PS71 values for Qx are very close to the "Maxwellian solution" . We assume that the motion 
of a is caused by Barnett relaxation. Consequently, we find perfect alignment of a with J 
for T s /T g = 0, with Qx decreasing monotonically to the Maxwellian solution as T s /T g — > 1. 
The behavior of Qx and Qj explains why our solutions for R are systematically larger than 
those of PS71.Q We conclude that the Barnett effect generally increases the efficiency of DG 
alignment. The enhancement is a factor of ~3 for T s /T g <C 1 but decreases rapidly as the 
dust temperature increases. For T s /T g ^ 0.3, there is little difference between our solutions 
and those of PS71. 

Figure 14 compares the 5 m dependence of our solutions with analogous solutions from 
PS 71, for the case of cold grains with T s /T g = 1/9. For small 5 m , our R values exceed those 
of PS71 by factors of 2-3. Our results appear to converge to those of PS71 for 5 m ^> 1. 
However, it is not clear whether the convergence occurs for other axis- and temperature 
ratios, there being no obvious physical explanation for this behavior. 



The essential difference between L97 and this paper is that the former replaces the 
random variables cos 2 9 and cos 2 (3 with their mean values in the Fokker-Planck equation. 
The correlation between 9 and /3, which was neglected in earlier studies, was introduced 
in L97 through varying the value of \J\. Therefore, it would be correct to write down the 
assumed dependence of R on Qj and Qx in L97 in the following way: 



where J|(cos 2 /3) expresses the conditional dependence of J on (3. We briefly discuss the 
perturbation procedure in Appendix C, as some points of it have been rectified after L97 
was published. In fact, this is the first actual use of the perturbation procedure developed 



To obtain R we start with the Maxwellian value of angular momentum (see LR97) 



5.2. 



L97 



R^QjQx(J\(cos 2 P)) 



(5-1) 



in L97. 




(5-2) 



12 This would not be the case if the "correlation factor", r, were different in the two studies. In fact, 
t = 0.6-0.7 for the PS71 data plotted in Fig. 13 (see PS71, Table 2), nearly identical to our mean r value. 
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corresponding to the temperature (see L95, L97): 

m T m + T s Oj_i 



(5-3) 



1 + Si-j. 



This enables us to find cos 2 (3 and cos 2 9 that are used to iterate the mean value of angular 
momentum and Qj and Qx- We find that the iterations quickly converge and that the mean 
value of grain angular momentum is much more sensitive to variations of (3 than of 9. 

The angles (3 and 9 are correlated within the perturbative approach through the value 
of the grain angular momentum. If we disregard small terms proportional to powers of small 
parameter 7, the value of the mean angular momentum that is used within the perturbative 
approximation is (see ( jUlUQ ) 



does (J 2 )i- But according to ( pl|) (J 2 )i determines the value of Qx and therefore cos 2 9 (see 



To calibrate the accuracy of the perturbative method, we have carried out a point-by- 
point comparison between the numerical R values compiled in Tables 2-6 and corresponding 
values obtained with the fifth-order iterative method (see Appendix C). We find that the 
mean error in the perturbative results is 16.8%, the sign indicating that the perturbative 
method (slightly) overestimates the true solution. Given the various idealizations in our 
grain model, this accuracy should be sufficient for many applications. 



Does our present paper mean that the theory of Davis-Greenstein alignment is complete? 
We see a couple of outstanding problems that we hope to handle in the near future: The 
present study describes oblate grains. Although there exists some evidence that aligned 
grains tend to be oblate (Hildebrand 1988) a study of prolate grain alignment is needed. 
Our solutions assume that the magnetic properties of grains are isotropic. This is an exellent 
approximation for paramagnetic grains, but ferromagnetic grains should be anisotropic (see 
DL98). This means that for sufficiently high 5 m we have to consider the anisotropyp| of the 
grain magnetic response. Therefore a quantitative study of anisotropy effects is necessary. 



13 In some instances, e.g., when a grain contains numerous superparamagnetic inclusions, the magnetic 
response can be isotropic despite high <5 m . This is probably an exception, however. 




eq. (H)). 



6. 



Future work 
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JS67 attracted the attention of researchers to the necessity of studying the magnetic 
response of grains. A recent study by DL98 clarified a number of issues, especially relevant to 
the high frequency magnetic response of candidate materials. We feel that the low frequency 
magnetic response deserves a similar study. We conclude that there is plenty of work to be 
done on Davis- Greenstein alignment. We will be lucky to finish this work by the half-century 
anniversary of DG51. 

7. Summary 

The principal results of this paper are as follows: 

1. We have calculated the efficiency of the Davis- Greenstein effect for oblate spheroidal 
grains of arbitrary axis ratio. Our description of the grain dynamics includes gas 
damping, magnetic damping and Barnett relaxation. 

2. We have calculated the alignment of ordinary paramagnetic grains in an "adiabatic 
approximation" that is accurate to O(e), where e ~ 10~ 3 is the ratio of the timescale 
for Barnett relaxation to the dynamical timescale for the angular momentum. 

3. For ordinary grains, we have tabulated the Rayleigh factor, R, on a grid of >1000 
combinations of the relevant parameters. We have estimated the accuracy of these 
results using various benchmark tests. The error in R is typically a few units in the 
third decimal place. 

4. The main difference between our paper and classical studies of DG alignment is the 
incorporation of Barnett relaxation in the present study. Including the Barnett effect 
enhances R by factors up to 3 for small dust-to-gas temperature ratios but makes 
little difference for T s /T g > 0.3. 

5. We have calibrated the accuracy of an approximate method proposed in L97. The 
mean error in values of R computed via the perturbative method is ~17%. 

6. We have derived an exact analytic solution for the alignment of superparamagnetic 
grains. The solution is based on an adiabatic approximation that is accurate to O (S^ 1 ). 

7. Our calculations on superparamagnetic grains show that DG alignment is only efficient 
if the grains are much colder than the gas. The commonly-held view that superparam- 
agnetism produces "perfect" alignment under realistic conditions is a misconception. 
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TABLE 1 
Dynamical Timescales 



Symbol Definition 

t ro t Rotation period 

tLar L armor period 

tear Barnett time 

t gas Gas damping time 

t mag Magnetic damping time 



Value (s) 

7 xlO-s pig T-^(a/b)VH^ 

2xl0 5 Ps , (x'^y 1 Bz\bl, 

4 x lO 6 o 2 K- 1 T- 1 (a/6) ( 1+aV62 ) 3 b 7 

4 X 1U p s0 A_ 13 J g l (1 _ a 2 /6 2) 0_5 

7 x 10 9 T^ 1 p Si0 nj 1 T"^ 2 (a/6) 6_ 5 
4xl0 12 p sfi KZl 3 BZ 2 5 b 2 _ 5 



Timescales for a homogeneous, oblate spheroid rotating about its symmetry axis with kinetic energy kT g /2 
in a gas of pure H2 (after RDGF93, Table 1). The quantity K = \" /lo depends on the uncertain grain 
composition. The scaling of K is appropriate for ordinary paramagnetic substances; for superparamagnetic 
grains, increase K by an uncertain factor <J 10 5 (Draine 1996). 
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TABLE 2 

Rayleigh Reduction factor for a/b = 0.1 



log 











.1 




.2 




.3 




.4 




.5 




.6 




.7 




.8 


-1 


.00 


0. 


.018 


0. 


.014 


0. 


.008 





.006 






















-0. 


.90 


0. 


022 


0. 


014 





,010 





.006 






















-0. 


.80 


0. 


030 


0. 


019 


0. 


.011 





.008 





.005 


















-0. 


,70 


0. 


035 


0. 


.023 


0. 


.015 





,009 





,006 


















-0 


.60 


0. 


.046 


0. 


029 





,018 





.011 


0. 


.007 


0. 


.005 














-0. 


.50 


0. 


.055 


0. 


035 





,021 





.014 





,009 


0. 


.006 














-0 


.40 


0. 


066 


0. 


040 


0. 


.025 





.016 


0. 


.011 





.007 


0. 


005 










-0. 


,30 


0. 


.083 


0. 


.049 


0. 


.029 





,019 


0. 


,012 





.009 


0. 


.006 










-0. 


.20 


0. 


.101 


0. 


.059 





,034 





.022 





.014 





.010 


0. 


.006 










-0. 


.10 


0. 


.119 


0. 


.068 





,039 





.025 





,016 





.011 


0. 


008 


0. 


005 











0. 


.142 


0. 


.078 


0. 


.045 





.027 


0. 


.019 





.013 


0. 


008 


0. 


.005 









,10 


0. 


169 


0. 


.090 


0. 


.050 





,032 





,021 





.013 


0. 


.009 


0. 


006 









.20 


0. 


.200 


0. 


101 





.056 





.035 


0. 


.022 





.015 


0. 


010 


0. 


007 






0. 


.30 


0. 


.228 


0. 


.112 





,061 





.038 


0. 


,025 





.017 


0. 


Oil 


0. 


007 






0. 


.40 


0. 


.263 


0. 


.121 





,066 





.041 





,026 





.017 


0. 


012 


0. 


.008 






0. 


.50 


0. 


300 


0. 


133 


0. 


.070 





.043 


0. 


.028 





.018 


0. 


013 


0. 


.008 


0. 


.005 





,60 


0. 


337 


0. 


140 


0. 


.074 





,045 





,029 


0. 


.020 


0. 


013 


0. 


.008 


0. 


,005 





.70 


0. 


.372 


0. 


149 





,078 





.047 


0. 


.030 





.020 


0. 


013 


0. 


009 





.005 


0. 


.80 


0. 


.411 


0. 


.154 





,081 


0. 


.050 


0. 


,032 





.021 


0. 


014 


0. 


009 





,005 


0. 


.90 


0. 


452 


0. 


.161 


0. 


.082 





.050 


0. 


.032 





.022 


0. 


.014 


0. 


.009 


0. 


.005 


1 


,00 


0. 


.490 


0. 


.167 


0. 


.086 





,051 


0. 


,033 


0. 


.022 


0. 


015 


0. 


.009 





,005 


1 


.10 


0. 


,527 


0. 


.170 





.087 





.052 


0. 


.034 





.022 


0. 


,015 


0. 


009 





.005 


1 


.20 


0. 


562 


0. 


.175 





,088 





.053 





,034 





.022 


0. 


015 


0. 


009 





,005 


1 


.30 


0. 


597 


0. 


.176 


0. 


.089 





.054 


0. 


.034 





.023 


0. 


015 


0. 


.009 


0. 


.005 


1 


.40 


0. 


629 


0. 


.179 


0. 


.090 





,054 


0. 


,034 


0. 


.023 


0. 


.015 


0. 


009 


0. 


.005 


1 


.50 


0. 


661 


0. 


.178 





,089 





.054 


0. 


.034 





.023 


0. 


014 


0. 


009 


0. 


.005 
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TABLE 3 

Rayleigh Reduction factor for a/b = 0.3 



log 











.1 




.2 




.3 




.4 




.5 




.6 




.7 




.8 


-1 


.00 


0. 


.015 


0. 


on 





.007 


























-0. 


.90 


0. 
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0. 


015 





,009 





.006 






















-0. 


.80 


0. 


026 


0. 


018 


0. 
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.007 


0. 


.005 


















-0. 


,70 


0. 


031 


0. 
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0. 


.013 





,008 





,006 


















-0 


.60 


0. 
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0. 
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.011 


0. 


.007 


















-0. 


.50 


0. 


.053 


0. 
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0. 


,019 


0. 


.012 





,008 


0. 


.005 














-0 


.40 


0. 
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0. 
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0. 
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.015 


0. 


.009 





.006 














-0. 


,30 


0. 


.079 


0. 


.046 


0. 


.027 





,018 


0. 


,012 





.008 


0. 


.005 










-0. 


.20 


0. 


.098 


0. 
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.013 





.009 


0. 


.006 
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.10 


0. 
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0. 
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,036 


0. 
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,015 


0. 


.010 


0. 
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0. 
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0. 


.074 


0. 


.042 





.026 


0. 
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.011 


0. 
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0. 


.005 
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0. 


.167 


0. 


.083 


0. 
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,028 


0. 
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.012 


0. 


.008 


0. 
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.20 


0. 
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0. 
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.014 


0. 


.009 


0. 
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0. 


.30 


0. 
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0. 
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,055 





.034 
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.015 


0. 
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0. 
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0. 


.40 


0. 


.260 


0. 
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,059 





.036 





.024 





.015 


0. 


010 


0. 


007 






0. 


.50 


0. 
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0. 
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0. 


.064 





.039 


0. 


.025 





.016 


0. 


011 


0. 


.007 









,60 


0. 
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0. 


.128 


0. 


.066 





,040 





.026 


0. 


.018 


0. 


012 


0. 


.007 









.70 


0. 


,375 


0. 
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,070 





.042 
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.018 


0. 


012 


0. 
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0. 


.80 


0. 


.412 


0. 


.141 





,073 


0. 


.043 


0. 


,028 


0. 


.019 


0. 


013 


0. 


007 





,005 


0. 


.90 


0. 


.451 


0. 


,148 


0. 
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.044 


0. 


.029 





.019 


0. 


013 


0. 


.008 


0. 


.005 


1 


,00 


0. 


.490 


0. 


.151 


0. 


.076 





,046 


0. 


,029 


0. 


.020 


0. 


013 


0. 
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,005 


1 


.10 


0. 


.526 


0. 
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.046 
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.020 


0. 


013 


0. 


008 


0. 


.005 


1 


.20 


0. 


562 


0. 
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.020 
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,005 
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0. 
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0. 
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0. 
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0. 
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013 
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.005 
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0. 
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.020 
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TABLE 4 

Rayleigh Reduction factor for a/b = 0.5 



Ts/T % 



o 


5 m 









.1 




.2 




.3 




.4 




.5 




.6 




.7 


-1 


.00 





.014 





.010 


0. 


.005 






















-0. 


.90 





.017 





,012 


0. 


.007 


0. 


005 


















-0. 


.80 





.022 





.015 


0. 


009 


0. 


005 


















-0. 


,70 





.032 


0. 


.020 


0. 
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0. 
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0. 


005 
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.60 





.039 
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013 
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.008 
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.006 
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,50 





.049 
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016 
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010 
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0. 


.005 
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.40 





.064 
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.035 


0. 


019 


0. 


012 


0. 


008 


0. 


005 










-0. 


,30 


0. 


.079 


0. 


.041 


0. 


.022 


0. 


014 


0. 


009 


0. 


006 










-0. 
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,048 


0. 


.026 


0. 


016 


0. 


.010 


0. 


007 





.005 






-0. 


.10 
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,056 


0. 


031 


0. 


018 


0. 


012 


0. 


.008 





.005 














.139 


0. 


.064 


0. 
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0. 


021 


0. 


013 


0. 


009 





.006 









.10 


0. 
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0. 


.072 


0. 
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0. 
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.041 
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0. 
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0. 
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0. 
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0. 
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0. 
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.008 


0. 


,005 


0. 
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.298 


0. 


.101 


0. 


051 


0. 


030 


0. 


019 


0. 


.013 





.009 


0. 


.006 





.60 





.335 


0. 


.105 


0. 


052 


0. 


.032 


0. 


.021 


0. 


.013 





,009 


0. 


006 


0. 


.70 





.371 





.111 


0. 


055 


0. 


032 


0. 


.020 





,014 





.009 


0. 


006 





.80 


0. 


.409 





,116 


0. 


.056 


0. 


034 


0. 


022 





,015 


0. 


.009 


0. 


006 





.90 





.453 





.120 


0. 


059 


0. 


034 


0. 


022 


0. 


.015 





.010 


0. 


.006 


1 


,00 





.492 





.124 


0. 


060 


0. 


035 


0. 


.022 


0. 


.015 





,010 


0. 


.006 


1 


.10 





.527 





.126 


0. 


061 


0. 


035 


0. 


,023 





.015 





.010 


0. 


006 


1 


,20 





.564 


0. 


,127 


0. 


.062 


0. 


036 


0. 


023 





.015 





.010 


0. 


006 


1. 


.30 





.598 





.129 


0. 


062 


0. 


.037 


0. 


023 


0. 


.015 





.010 


0. 


.006 


1 


.40 


0. 


.631 


0. 


,131 


0. 


.062 


0. 


036 


0. 


023 


0. 


.016 





,010 


0. 


006 


1. 


.50 





.663 





.131 


0. 


063 


0. 


037 


0. 


.024 


0. 


.015 





.010 


0. 


006 



-25 - 



TABLE 5 

Rayleigh Reduction factor for a/b = 0.7 



log 


5 m 









.1 




.2 




.3 




.4 




.5 




.6 


-1 


.00 


0. 


013 


0. 


.007 






















-0. 


.90 


0. 


020 





,008 





.005 


















-0 


.80 


0. 


027 


0. 


.012 





.006 


















-0. 


.70 


0. 


.032 


0. 


.015 





,007 


















-0 


.60 


0. 


.042 





,018 





.009 





.005 














-0. 


.50 


0. 


.053 





,020 





.011 





,006 














-0 


.40 


0. 


064 


0. 


.025 





.012 


0. 


.007 


0. 


.005 










-0 


.30 


0. 


.080 


0. 


.030 





,015 


0. 


,008 





.005 










-0 


.20 


0. 


096 





.035 





.017 





.010 





.006 










-0. 


.10 


0. 


.118 





,039 





.020 





,012 


0. 


.007 





,005 











0. 


.141 


0. 


.046 





.021 





.012 





.008 





.005 






0. 


.10 


0. 


.167 


0. 


.050 





,024 





,014 


0. 


.009 


0. 


,006 









.20 


0. 


195 





.056 





.027 


0. 


.015 





.010 





.007 


0. 


005 





.30 


0. 


.227 





,063 


0. 


.029 





,017 





.011 





,007 


0. 


005 





.40 


0. 


261 





,066 





.030 


0. 


,018 





.011 





,008 


0. 


005 





.50 


0. 


296 


0. 


.071 





.032 


0. 


.019 





.012 


0. 


.008 


0. 


.005 


0. 


.60 


0. 


335 


0. 


.073 





,033 





,020 


0. 


.013 


0. 


,008 


0. 


.006 





.70 


0. 


376 





,077 





.035 


0. 


.020 





.013 





.009 


0. 


006 





.80 


0. 


.413 





,079 





.036 


0. 


,021 





.013 





,009 


0. 


006 





,90 


0. 


.451 


0. 


.082 





.037 


0. 


.021 





.014 


0. 


.009 


0. 


.006 


1 


.00 


0. 


.490 


0. 


.083 





,038 





,022 





.014 


0. 


.009 


0. 


006 


1 


.10 


0. 


.528 





.084 





.039 





.022 





.014 





.009 


0. 


006 


1 


.20 


0. 


563 





,085 





.039 





,022 





.015 





,009 


0. 


006 


1 


.30 


0. 


.598 


0. 


.087 





.039 


0. 


.023 





.015 


0. 


.010 


0. 


.006 


1 


.40 


0. 


.630 


0. 


.088 





,040 





,023 





.015 


0. 


,010 


0. 


006 


1 


.50 


0. 


663 





,088 





.040 





.023 


0. 


.015 


0. 


.010 


0. 


006 
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TABLE 6 

Rayleigh reduction factor for a/b = 0.9 



log<5 m 





.1 


.2 


.3 


.4 


-1.00 


0.015 














-0.90 


0.021 














-0.80 


0.026 














-0.70 


0.035 


0.005 











-0.60 


0.042 


0.007 











-0.50 


0.056 


0.008 








-0.40 


0.068 


0.009 








-0.30 


0.082 


0.011 


0.005 






-0.20 


0.100 


0.013 


0.006 






-0.10 


0.119 


0.015 


0.006 









0.141 


0.016 


0.007 






0.10 


0.169 


0.018 


0.008 


0.005 




0.20 


0.198 


0.020 


0.009 


0.005 




0.30 


0.227 


0.022 


0.010 


0.005 




0.40 


0.261 


0.023 


0.010 


0.006 




0.50 


0.300 


0.024 


0.011 


0.006 




0.60 


0.335 


0.025 


0.011 


0.006 




0.70 


0.372 


0.026 


0.012 


0.007 




0.80 


0.414 


0.027 


0.012 


0.007 




0.90 


0.450 


0.028 


0.012 


0.007 


0.005 


1.00 


0.489 


0.029 


0.013 


0.007 


0.005 


1.10 


0.526 


0.029 


0.013 


0.007 


0.005 


1.20 


0.563 


0.030 


0.013 


0.007 


0.005 


1.30 


0.598 


0.030 


0.013 


0.007 


0.005 


1.40 


0.631 


0.030 


0.013 


0.008 


0.005 


1.50 


0.661 


0.030 


0.013 


0.008 


0.005 
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TABLE 7 

Analytic Solution for Superparamagnetic Grains 
vs. Numerical Solution for 6 m = 100 



a/b 




Analytic 


Numerical 


0.1 


0.1 


0.1887 


0.1857 


0.1 


0.2 


0.0951 


0.0940 


0.1 


0.3 


0.0571 


0.0559 


0.1 


0.4 


0.0370 


0.0366 


0.5 


0.1 


0.1369 


0.1353 


0.5 


0.3 


0.0388 


0.0380 


0.5 


0.5 


0.0166 


0.0165 



TABLE 8 

Rayleigh reduction factor for superparamagnetic grains 



a/b 










T s /T g 










.1 


.2 


.3 


.4 


.5 


.6 


.7 


.8 


.9 


0.1 


.1887 


.0951 


.0571 


.0370 


.0247 


.0164 


.0106 


.0061 


.0027 


0.3 


.1688 


.0839 


.0502 


.0324 


.0217 


.0144 


.0093 


.0054 


.0024 


0.5 


.1369 


.0657 


.0388 


.0250 


.0166 


.0111 


.0071 


.0041 


.0018 


0.7 


.0916 


.0417 


.0242 


.0155 


.0103 


.0068 


.0044 


.0025 


.0011 


0.9 


.0316 


.0138 


.0080 


.0051 


.0034 


.0023 


.0014 


.0008 


.0004 
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A. Diffusion coefficients in equation (3-1) 

A.l. Gas damping 

The diffusion coefficients for gas damping were derived elsewhere (LR97; L97). The 
(dimensionless) mean torque has body frame components 



h (tjtA jI 



A lv = -h (r±/r„) J* 



and 



A b — — V 



The diffusion tensor is diagonal in the body frame with 

B h s , xx = (rjL/r,,) (i + T s /T g ) 



R — R 

g,yy ~ 



and 



nb 



(i + r./r g ). 

The dimensionless shape factors 

rn(e) = ^ { 3 + 4(1 - e 2 Me) - e" 2 [l - (1 - e 2 fg(e)} } 

and 



r ± (e) = -{7-e 2 + (l-e 2 )Ve) + (l-2e 2 



1 + e 



are weak functions of the eccentricity, e = J 1 — a 2 /b 2 , with 

qf(e) = — In ( — — — 
yK J 2e VI -e 



(Al) 
(A2) 

(A3) 

(A4) 
(A5) 

(A6) 
(A7) 



l-(l-e 2 ) 2 s(e)]]}, (A 



(A9) 



The diffusion coefficients in equation ( |3-1| ) are evaluated in the inertial frame. To obtain 
these coefficients from expressions ([A1|) -(|AT|) we (i) transform components from the body 
frame to the inertial frame; and (ii) average the resulting expressions over the rapidly- 
changing angles and ip (see L97). The calculations are straightforward but tedious and 
here we simply give the results. The mean torque has inertial frame components 



A g:i = — ( 1 + A sin 2 Ji i = x,y,z, 



(AlO) 
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where 



a = h (r ± /r,|) - 1. 

The diffusion tensor has inertial frame components 



5, 



g,xx 



Ul+ T s /T g ) [ (l + cos 2 p) P 1 + sin 2 (3 P 2 



Bg,yy — B gyXX) 



and 

where 

and 

are functions of 9 and 



B &zz = (1 + T s /T g ) [ sin 2 /3 Px + cos 2 /3 P 2 
Pi = 1 - 7 (l - sin 2 0/2) 
P 2 = (l - 7 sin 2 0) 

7 = i-r ± /r,|. 



(AH) 

(A12) 
(A13) 

(A14) 

(A15) 

(A16) 

(A17) 



A. 2. Magnetic damping 

The mean torque was derived by DG51 (see DG51, eq. [81]). In dimensionless units, 
the inertial frame components are 



^m,i — Z( y 9)5 ni J x , 

Am,y = ~Z (9) 5 m J y , 



and 
where 



0, 



Z(0)= 1 + (h - 1) sin 2 9 



(A18) 
(A19) 

(A20) 

(A21) 



The diffusion tensor for paramagnetic relaxation in spheroidal grains does not appear 
in previous papers on the DG effect, which either ignored thermal fluctuations (e.g., DG51) 
or used the diffusion tensor for spheres (JS67). However the diffusion tensor is uniquely 
determined by the principle of detailed balance, which requires the probability current at 
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each point in phase space to vanish in thermodynamic equilibrium. Setting the current along 
the x direction to zero yields 

i d 

^m.z /* ~ 2~dJ~ ^ m,xx ^' = ^' (A22) 

where 

r = c ^{^) (A23) 

is the thermal equilibrium distribution of the dimensionless angular momentum and C is a 
normalization constant. Consistent with the symmetry of the problem, we have assumed 
that the off-diagonal components of the diffusion tensor are zero. The unique solution of 
eq. (|A23|) with B m xx finite for all J x is 

B m , xx = 2 (T s /T g ) 5 m . (A24) 

The analogous conditions for the y and z directions yield 

Bm,yy B mxx (A25) 

and 

B m , zz = 0. (A26) 
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B. Analytic solution for superparamagnetic grains 



B.l. Adiabatic elimination 

We seek an approximate solution to equation ( |3-8D of the form 

fext {^Jxi Jy, J z) fxy (Jxy Jy\Jz) fz \'^z) 

subject to the normalization conditions 

/+oo r+oo 
/ fxy \ Jxi Jy\Jz) dJ x dJy 1 

-co J —oo 

and 



+ CO 



fz{J z ) dJ z 



(Bl) 



(B2) 



(B3) 



We require our approximation to be accurate to O This means that we may neglect the 

terms of order unity in comparison with the terms of order # m in expressions ( |3-12|) - (|3-14 ) 
and ( |3-19|) -( j3-21| ) for the diffusion coefficients. For the linear coefficients, this gives 



A x » - [1 + (h - l)a] J x 5 m = A x 6 m , 
A y « - [1 + (h - l)a] J y 6 m = A y 6 m , 



and 



A z = -(1 + Act) J z . 
The analogous approximations for the quadratic diffusion coefficients are 

B xx w 2 (T s /T g ) 5 m = B xx 8 m , 

B„ 

and 

By-, 



'yy 



Bxx "m; 



(B4) 
(B5) 

(B6) 

(B7) 
(B8) 

(B9) 



1+T s /T g ) [sin 2 /5P 1 (J) + cos 2 /3P 2 (J) 

where -Pi (J) and P<z{J) are defined in eqs. (|A15|) and ( |A16| ). The physical content of our 
approximation is that the torques perpendicular to B are provided entirely by magnetic 
dissipation; evidently this is accurate to order t mag /t gas . 

and (|B~7|)-(|B~9|) into equation 

(BIO) 



To proceed, we substitute expressions ( |B"T 
(|3-£j|). After a little algebra, the latter becomes 



fz^xyfxy ~t" ^ m ^zfxyfz 0) 
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where 

c) ~ d ~ 1 d 2 ~ Id 2 



and 



dJ x A ' X ' dJy^^ 2dJ 2 x BxX ' + 2dJ 2 y BxX ' (BU) 

d - Id 2 - 

c - s -w. A -- + 2Bjf B -- (B12) 



B.2. Solution for f xy 

Neglecting the term of order 5" 1 in equation ( piO| ) gives 

£xyfxy = 0, (B13) 

which is a Fokker-Planck equation for f xy . After transforming^ to polar coordinates 

1/2 



and 



(p = tan -1 (J y /J x ) , 



equation ( |B13j ) becomes 




dJ\ 



A±f± — -B xa 



2 xx dJi 



d 



~2 \ B xx/J± 



where 
and 



f± = J±f, 



A± = -[l + (h - l)a] J ± + 



B , 



2Ji 



(B14) 
(B15) 

(B16) 

(B17) 
(B18) 



Equation ( [B16D will be satisfied if we assume that f± is independent of 0, consistent with 
the symmetry of the problem, and require 



A±f± — -Bxx-T-r- : ~ 0. 



2 xx dJ 



(B19) 



Equation (|B19|) states that there is no probability current perpendicular to the B direction 
in equilibrium ( "zero- current solution" for 



1 The transformation is straightforward but tedious and only the relevant results are given here. For 
details, see any text on statistical mechanics (e.g., Risken 1984). 
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Solving eq. ( |B19|) for f± and using eq. (|B17| ), we find 

f xy = C x exp (-$) , (B20) 
where C x is determined by eq. ( |B"2"D and 

$ (J x , J y \J z ) = Q>j J q J± [l + (h-l)a (J')] dJ' ± . (B21) 

Since does not involve differentiation with respect to J z , the latter appears^ only as a 
parameter in the solution for f xy . In physical terms, this means that the distribution of J x 
and J y relaxes rapidly ( "adiabatically" ) to the equilibrium distribution that would obtain if 
J z were frozen at its instantaneous value. Of course, this is appropriate when t mag <C t gas . 



B.3. Solution for f z 

Given the solution for f xy from the preceding section, one can can find f z by setting the 
term O in equation ( [B10|) to zero. When eq. (|B12|) is used to expand C z in terms of 
derivatives, this statement becomes 



d 
d7, 



A-zfxyfz 



+ 



1 d 2 



2d J 2 



B zzf xyf z 



0. 



(B22) 



If the last equation is multiplied by dJ x d.J y and the resulting expression is integrated over 



all J x and J y , the result is 



where 



and 



d 
dJ z 

A z 
B 7 , 



A z fz 2 (^Bzzfz^j 



r+co r+oo a r A 1 A 1 
J-oo J—oo ^-zjxy mux uui 



(B23) 



(B24) 



J^J^Bzzf xy dJ x dJ y . (B25) 

Equations (P23|) -( [B25|) say that f z satisfies a one-variable Fokker-Planck equation wherein 
the diffusion coefficients A and B have been averaged over J x and J y in the obvious way. 
The preceding discussion shows that this procedure is accurate to O (S^ 1 ). 



15 Recall that a depends on J 2 , hence J z . 
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The zero-current solution for f z is elementary: 

fz = C z exp 

where 



This completes the solution. 



J - f-2A, 



o V B ? 



dJ' 



(B26) 
(B27) 



C. Perturbative scheme 

It was shown in L97 (see also Errata in Lazarian 1998) that calculations of (Qj) can be 
performed perturbatively with 



(Q^^cos 2 /?,)-! 



where 



(cos 2 /3 m ) 



K 2 
4 \+i 



T 

n< = 4^ + 



1 - K 2 

7 



=arcsm 



i+l 



-QxiQ 



(Cl) 



(C2) 



(C3) 



T av is given by ( p-3| ) and T m is the mean rotational temperature of grain in the absence of 
paramagnetic relaxation. For the sake of simplicity this temperature is frequently chosen to 
be equal 0.5 (T s + T g ), but such a choice is uncertain (see Purcell 1979). 



The internal alignment measure Qxi can be found as (LR97, L97) 

3 (ft cos 2 sin Of (6) d0 l\ 



Q 



Xi 



and 



2 V /o* sin 0/(0) d0 3 J 
/(0) = const x sin exp [-Efi(0)/kT t 

(J 2 ) 



2I Z 



l + (h- l)sin 2 



where (J 2 )i is to be found using the distribution function for the i-th iteration is 

/ = fx fy fz ~ const x exp (-a J 2 ) 

where 



(C4) 
(C5) 

(C6) 
(C7) 
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;i + sin 2 ^[/ir ± /rn - 1]) j 2 (i - cos 2 a [i - k 2 ]) v2(/3 i ,e i )T m 

Ot— ; : 1— — ; , 



1 



1 



1 + cos 2 0) ( 1 + cos 2 0) + sin 2 /3 sin 2 



(C9) 



In the perturbative approach all the values obtained at the i stage are treated as con- 
stants and therefore integration provides 



poo I 

( J)i = const x J 2 exp(aJ)JdJ = — 
Jo 2a 



(CIO) 
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FIGURE CAPTIONS 

Fig. 1 - - Definition of the body frame, ja; 6 , y b , z h \, with z b parallel to the symmetry 
axis, a, and the remaining basis vectors oriented arbitrarily in the equatorial plane of the 
grain surface (the lightly-shaded disk). The orientation of a with respect to the angular 
momentum, J, is specified by angles 9 and ip. The "J frame" basis, [x J ,y J , z J }, is used 
only to carry out transformations between the body and inertial frames (see Fig. 2). 

Fig. 2 — Definition of the inertial frame, {x, y, z}, with z parallel to the interstellar mag- 
netic field, B, and the other basis vectors oriented arbitrarily in the plane perpendicular to 
B. Angles j3 and <fi specify the orientation of J with respect to B. 



Fig. 3 — The cross sections for polarized light are defined with respect to the observer 
frame basis [see eqs. ( |2-5| ) and ( |2-6| )1, with z° the direction of propagation and y° along the 
projection of B onto the plane of the sky. 

Fig. 4 — In the adiabatic approximation, statistics of the 9 distribution depend only on £ 2 [a 
dimensionless function of J 2 , see eq. (|3-4|) ], the grain shape and the dust-to-gas temperature 



ratio. Shown here are the mean values of cos 2 9 (short dash), sin 2 9 (long dash), and qx 
(solid). 

Fig. 5 — The angular momentum diffusion coefficients, and hence the Rayleigh reduction 
factor, depend on the grain shape. For a homogeneous grain, the dependence is only via the 
functions of a/6 shown here: A (short dash), 7 (long dash), and h — 1 (solid). 

Fig. 6 — Histogram: distribution of the errors in Qj, computed by comparing the results of 
100 identical numerical trials with exact solution for spherical grains (see §4.1). The mean 
and rms errors in Qj are —1.5 x 10~ 4 and 1.2 x 10 -3 , respectively. Dashed curve: Gaussian 
distribution with the same mean and standard deviation as the histogram. 

Fig. 7 - - Comparison between Qj values computed numerically (open circles) and the 
predictions of the perturbative approach described in Appendix C. The calculations shown 
here are for the case T s /T g = 0, a/b = 2/3. 

Fig. 8 — Similar to Fig. 6 but showing the errors in Qx for a special case, T s /T g = 1, where 
an exact solution for Qx exists (see §4.1). The mean and rms errors in Qx are —5 x 10~ 4 
and 5 x 10~ 4 , respectively. 

Fig. 9 — (a) Solid curves: values of R (top), Qj (middle) and Qx (bottom) determined 
by numerical integration of the Langevin equations, plotted vs. the axis ratio. Values of 
the dust-to-gas temperature ratio and magnetic damping parameter are indicated. Dashed 
curves: product of QxQj from the numerical calculations (top), the exact value of Qj for 
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a spherical grain with the same T s /T g and S m (middle) and the exact value of Qx for a 
Maxwellian angular momentum distribution in the grain frame (bottom), (b) Similar to (a) 
but for 5 m = 10. 

Fig. 10 Similar to Figure 9 but showing the dependence on the magnetic damping 
parameter with a/b and T s /T g fixed. 

Fig. 11 — Similar to Figure 9 but showing the dependence on the dust-to-gas temperature 
ratio with a/b and S m fixed. 

Fig. 12 — Dependence of the correlation function [see eq. Q4-6Q 1 on the magnetic damp- 
ing parameter. Solid curve: numerical calculations with the other parameters fixed at the 
indicated values. Dashed curve: approximation to p obtained by replacing the angular mo- 
mentum distribution with the angular momentum distribution for spheres with the same 
T s /T g and S m (see text). 

Fig. 13 — Comparison between the results of PS71 (filled circles) and our numerical results 
(solid). The dotted curves depict Qj for spherical grains with the same S m and T s /T g (middle 
panel) and Qx for a Maxwellian angular momentum distribution in the body frame (bottom 
panel). 

Fig. 14 — Similar to Fig. 13 but showing the dependence on 5 m . 
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